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We present a phenomenological interaction with a scale factor power law form which leads to 
the appearance of two kinds of perturbed terms, a scale factor spatial variation along with per¬ 
turbed Hubble expansion rate. We study both the background and the perturbation evolution 
within the parametrized post-Friedmann scheme, obtaining that the exchange of energy-momentum 
can flow from dark energy to dark matter in order to keep dark energy and dark matter den¬ 
sities well defined at all times. We combine several measures of the cosmic microwave back¬ 
ground (WMAP9-bPlanck) data, baryon acoustic oscillation measurements, redshift-space distor¬ 
tion data, JLA sample of supernovae, and Hubble constant for constraining the coupling con¬ 
stant and the exponent provided both parametrized the interaction itself. The joint analysis of 
Planck -I- WMAP9 -I- BAO -l-RSD -I- JLA + HST data seems to favor large interacting coupling, 

= 0.34403427lQ )^ggQy 353 at 1 ct level, and prefers a power law interaction with a negative ex¬ 
ponent, thus jS = —0.50863232lQ’4Qg23857 at 1 cr level. The CMB temperature power spectrum 
indicates that a large coupling constant produces a shift of the acoustic peaks and affects their 
amplitudes at lower multipoles. In addition, a larger /3 exponent generates a shift of the acoustic 
peaks, pointing a clear deviation with respect to the concordance model. The matter power spec¬ 
trum are sensitive to the variation of the coupling constant and the /3 exponent. In this context, 
the interaction alters the scale of matter and radiation equality and pushes it away from the present 
era, which in turn generates a shift of the turnover point toward to smaller scale. 

PACS numbers: 98.80.-k, 98.80.Es 


I. INTRODUCTION 

Our current view of the Universe is based on the large 
amount of cosmological observations provided by several 
surveys m. 0 : 0 - In particular, the first released data 
by Planck mission seem to indicate that our Universe is 
filled with two dark components representing the 95 % of 
the cosmic budget 01 : 0 : 0 : in fact, these components 
are non-baryonic in nature. On one side, dark energy 
can be treated as a fluid with negative pressure which vi¬ 
olates the strong energy condition, counter-balancing the 
net gravity effect on the large scale by driven the Universe 
toward an accelerated expansion era, leaving behind an 
earlier non-accelerated epoch. However, it was not dis¬ 
covered yet what is the microscopic mechanism which 
triggers the expansion of the Universe and therefore there 
is not a single fundamental theory which shows how fast 
this transition takes place or when exactly occurs. On 
the other side, dark matter can be accommodated as a 
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pressureless fluid which helps to cluster matter and af¬ 
fect in this way the whole distribution of galaxies in the 
Universe 0: 0: 0- Such component leaves its imprint 
in the power spectrum of clustered matter or through 
gravitational leasing of mass distribution 0], 0. A fur¬ 
ther exploration of the clustered matter dynamic needs 
to include observations that account for the behavior of 
dark matter perturbation and its clustering effect on sub- 
Hubble scale. A useful manner to achieve this goal is 
by redshift space distortion (RSD) measurements of fas 
quantity at different redshifts, where / is the growth rate 
factor and erg is rms density contrast within an 8 Mpeh”^ 
sphere, while h stands for the dimensionless Hubble pa¬ 
rameter m, 01]: 01I- Given RSD measurements are 
related with peculiar velocities of galaxies, the use fas 
data is considered as a promising way to test a cosmo¬ 
logical model beyond the background level, helping to 
break the degeneracy between different parameters 03], 
03], 05], 06]. 

A pressing issue of modern cosmology relies on the ex¬ 
istence of non vanishing interaction between dark matter 
and dark energy. This phenomenon is not only consis¬ 
tent with recent Planck cosmological data 07], 08], 09], 
but also seems to be possible at theoretical level when 
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coupled scalar fields are considered [20], jH], |12|, HO). 
A signature of the exchange of energy in the dark sec¬ 
tor can be searched in different manners, for instance, 
a starting point is to explore the behavior of dark en¬ 
ergy density at early time |24j along with the amount 
of dark energy at the recombination epoch |5], [25] . 
[26| . Another interesting route is to examine the im¬ 
pact of the cosmological constraints over an interacting 
dark energy model when a dynamical probe such as RSD 
measurements are included [19j . As an exmple, a geo¬ 
metric (background) test based on CMB -|- BAO -I- SNe 
data for an interaction proportional to dark matter den¬ 
sity seems to avoid large coupling, leading to a cou¬ 
pling constant = 0.0013 ± 0.0008 at Ict level [27] . 
The inclusion of RSD measurements improves a little 
bit this constraint provided the joint statistical analysis 
performed with CMB -|- BAO -f SNe -I- RSD data leads to 
= 0.00140j'jQ;gQggg at la level [28]. Even though the 
latter result differs from the former one, we notice that 
the disagreement is not bigger than 0.022%. In the light 
of previous outcome, one might wonder if the small value 
of the coupling constant can be extrapolated as a gen¬ 
eral output of interacting dark energy models when dark 
matter and dark energy are both treated as fluids. Con¬ 
trary to the previous results, Salvatelli et al. found that a 
piecewise vacuum interaction, with a large negative cou¬ 
pling constant, is favored at late time by the joint analysis 
of Planck -|- Union2.1 -|- RSD -|- BAO data [23]■ The lat¬ 
ter result would show that the coupling constant strenght 
depends on several factors, that is, it is likely linked with 
the specific form of the interaction, the data selected for 
implementing a statistical analysis, and above all it is not 
only connected with the underlying framework used for 
treating the dark components [28] . 

In this work, we are going to explore an alternative in¬ 
teraction between dark matter and dark energy which is 
parametrized as a product of the Hubble function and a 
scale factor power law up to a coupling constant. We will 
include linear perturbations into the cosmic constraint 
with the help of the parametrized post-Friedmann (PPF) 
formalism EO], m. m for avoiding non-adiabatic insta¬ 
bilities which could happen within the framework of an 
interacting dark sector [27], [33], [31], [33]. We are go¬ 
ing to assess how the interaction affects the background 
evolution along with its impact on the perturbation equa¬ 
tions in the IPPF scheme [27], m- We will analyze how 
the cosmic parameters are affected by this exchange of 
energy taking into account that the perturbed fluid equa¬ 
tions will admit a perturbed interaction which have both 
spatial variation of the scalar factor along with perturbed 
Hubble expansion rate. We will consider a momentum 
transfer potential corresponding to an interaction vector 
proportional to dark matter velocity in order to avoid 
a violation of the weak equivalence principle. We will 
perform a MCMC statistical analysis by modifying Cos- 
moMC package [36]. Our analysis will rely on different 
observational probes such supernovae la from JLA sam¬ 
ple, BAO, HST, Planck -I- WMAP plus a dynamical test 


based on RSD measurements. We also will explore the 
TT power spectrum of cosmic microwave radiation and 
the power spectrum of clustered matter by running a 
modified version of the CAMB code m, focusing on the 
behavior of both magnitudes when the coupling constant 
varies. 


II. BACKGROUND AND PERTURBATION 
EQUATIONS 

A. Background 

We assume a spatially flat Friedmann-Robertson- 
Walker background described by the metric = 
a^diag[—1, 7 ij] with local coordinates where the 

conformal time is defined as drj = a~^dt, and stands 
for the the spatial metric. Einstein equations, = 
with K = 8nG, lead to Friedmann constraint: 

— ~^{PT + Px) ■ (1) 

(X o 

As is customary, the continuity equations for total mat¬ 
ter (dark matter, baryons, neutrinos) and effective dark 
energy are given by 

p't = -3(pt +Pt) + ^ , 

p'x = —^{Px +Px) + ^ , (2) 

when dark matter and effective dark energy are coupled. 
For simplicity, we assume a phenomenological interaction 
that can be parametrized in terms of two constants only 
and that the form of the coupling between dark compo¬ 
nents is of a power law form for concreteness: Qx = HQj; 
with Qx = —^cPcoa^- We solve background equations for 
the interacting dark sector [38] : 


Px — PxQ^ 


— 3(l+lHa:) 


^cpcO 


3(1 + Wx) + P 


Pc — PcO^ 


— 3(1+iHc) 


^cpc{) 


3(1 + Wc) ~\~ P 


a^. 


(3) 


As an illustration, we write the cold dark matter density 
as the product of its density by its mass, Pc = ncfricia), 
and we realize that the mass of cold dark matter particle 
is no longer constant but it varies with the cosmic time: 

TOc (a) = pco (1 + || ■ (4) 


At this point we would like to emphasize that dark mat¬ 
ter and dark energy are both treated as fluids because 
the true nature of these dark components is not known 
yet provided the degree of freedoms which describe them 
have not been discovered so far. Furthermore, an interac¬ 
tion like we have proposed might not have a simple iden¬ 
tification in terms of a coupled quintessence model [2D], 
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m, m, m- Now, we focus on the possible choices 
that the parameter space would admit and make some 
comments about these cases. In the case of Wc = 0 and 
> 0, the parameter space is defined through the con¬ 
ditions /3 > —3 and Wx < —(1 + /3/3) < 0, which allow 
us to have positive energy density at all times. Another 
possible branch corresponds to taking Wc = 0 and < 0 
in the background equations ([^, however, this case leads 
to an unphysical conditions Wx > 0. In short, the in¬ 
teraction proposed here only allows to have an exchange 
of energy-momentum from dark energy to dark matter 
(^c > 0) in order to keep the energy densities as posi¬ 
tive definite quantities for all times. Besides, the vanilla 
model can be obtained when the coupling constant van¬ 
ishes together with the choices Wc = 0 and Wx = —1. 


B. Perturbation 

As we are interested in analyzing the behavior of scalar 
linear perturbation, we have to perturb the Einstein 
equation and balance equations around a FRW back¬ 
ground. To do that, we conveniently write the perturbed 
variables as a linear combination of the eigenfunction 
(plane-wave) Y (x) of the Laplace operator [39] , |4^ . At 
the end, the perturbed metric is characterized by four 
functions called potential A, shift B, curvature H^, and 
shear Ht'- 


Sgoo = -a^{2AY ), Sgo, = -a^BY ,, 

5g,, = a\2HLY-f,j + 2 HtY,, ). (5) 


momentum: 

= Q’^j, ^Q-, = 0. (8) 

I—x,c 

Then, we have a perturbed continuity and a perturbed 
Navier-Stokes equation for each fluid 


+ Api) -|- {pi Pi){k}iVi + 3iL^) 


_ AQi — ^Qi 
~ H ’ 

[a^{pi + Pi){Vi - B)]' 
a'^kn 


Ap* -t- -CKPiA^ 


= lmv-VT)+f.]. 


{pi -|-pi)A 

(9) 


For scalar perturbations ([^, the perturbed velocity 
vector becomes u,^ j = a[ — {1 + AY); [Vi — B)Yj\ , which 
means the four vector interaction can be written as 


Q.i = a{- Qi{l + YA) - SQiY; [Fi + Qj{V - B)]Yi) , 

( 10 ) 

where 6Qi and Fj stand for the energy transfer pertur¬ 
bation and the intrinsic momentum transfer potential of 
/-fluid, respectively. The conservation of total energy- 
momentum for dark components implies that 

Qc — Qx I Bq — Fx , 5Qc — ^Qx' 


C. The PPF method 


Here Yi stands for the covariant derivative of the eigen¬ 
function T(x). The perturbed energy-momentum tensor 
involves isotropic and anisotropic pressure perturbations, 
density, and velocity perturbations: 

ST°o = -Sp, 6To^ = -{p + p)vY\ 

ST) = SpYS) + pUY] . (6) 

The Einstein equations can be written as 


« i ‘ ff i ® 


V — B 


2H^CKk^ 

Sp + 3{p+p) 

kn 

•> 

A K ( 

B 

H), 


3 (ai/)2 V 

kn 

k^ 

K, . 

. V — B 



+ W k„ ■ 




(7) 


where the prime refers to a logarithmic derivative, ku = 
{k/aF[), ck = 1 —SAT/A:^, and K is the spatial curvature. 
In order to further illustrate how the perturbation theory 
works, we need to propose a covariant four-vector inter¬ 
action, Q'^j, which encodes the transfer of energy and 


To apply the PPF method we must make some as¬ 
sumptions about the behavior dark energy perturbations 
on super-horizon scales, that is, the effective dark energy 
contribution in the large scale limit ikn —t 0) must be 
accommodated in terms of a single function /^(a) |5T] : 

tSi , (11) 

Here we implicitly assumed to be working in the co¬ 
moving gauge {B = Vt, FIt = 0) where we called A = ^, 
Hi^ = (. Further, the derivative of curvature perturba¬ 
tion becomes 


,. p[Qc{V — Vt) + fc] + ^Pt — ^ckPtAt 

hm 4 = — -—^-- 

[Pt + Pt) 

— ^^hVt + -CKfcknVT , (12) 

so the PPF method deals with ~ 0 in the co-moving 
gauge, consequently dark energy fluctuations become im¬ 
portant at second order in the co-moving wave number 
ku (cf.[3T| ) provided knYr = 0{k‘j^(). From now on, 
we are going to neglect any contribution from anisotropic 
pressure terms for total matter and effective dark energy, 
Ht = Hj; = 0. 
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The behavior of the gravitational potential is related 
with the evolution of curvature perturbation by means 
of $ = C + VtI^h- In fact, the potential must fulfill a 
Poisson-like equation in the quasi-static limit: ^ 

1) = Ko?/S. tPt/‘^ k'^CK ■ The PPF prescription relies on 
the simple fact that the behavior of curvature perturba¬ 
tions on the two different scales must be linked by means 
of a single degree of freedom encoded in the P function 
m as follows: 


$-tr = 


2cKk^ 


^tPt ■ 


(13) 


In the large-scale limit, we obtain P' = S' — P, where the 
source term can be recast as 


Ka^ / 3a 


1 


5 =— — [Q,(l/ - Vt) + F,\ + — [AQ, - iQ,\ 


2k^ Hck 

+ VVfcff [-/c(Pt + Pt) + {Px + Px 


(14) 


In order to link the behavior of perturbed quantities at 
different scales, we must introduce a master equation for 
P: 


(i + c^fc|^)[p' + r + c^fc^r] = s, (15) 

where the condition cr^ = "H determines the transition 
between large scale regime and the quasi-static phase, 
that is, the PPP method introduces a new parameter 
called cr- In regard to the initial condition for P, we take 
a vanishing P at small scale factor provided the source 
term goes to zero in this limit. 

So far we have been working with the function P in 
order to connect two separated scales, but we have not 
said anything about how we deal with dark energy per¬ 
turbations in the previous method. The crucial fact of 
the PPF method is that forces dark energy density per¬ 
turbations to remain smaller than the dark matter per¬ 
turbation m- In doing so, one does not give a closure 
relation between dark energy pressure perturbation and 
dark energy density perturbation because one would like 
to avoid the appearance of large-scale instabilities m, 
pS] . As a result, the perturbed dark energy density and 
peculiar velocity are both derived once the new dynami¬ 
cal function is known (cf. [21], jUj): 

Px^x ^\Px “t” Px) 7 — 2 ^ 

K'a^Vx{px + Px) = - P - P'] 

F[a) V 

2 

+ ^/c(PT +PT)PTfcff) + na^Vripx+Px ), (17) 

where F(a) = 1 -I- ^ 2 k^^iPT +Pt)- Here we are keeping 
only the leading terms in the source expression for small 
comoving wave number kn, and we are calling this result 

So [Hj- 


In this study, we will be concerned with the physical 
consequences coming from the IPPF method when an 
interaction four-vector can be written as = HQ^u^ 
with Qx = —and Qx = aQ° = HQx- In this 
case, the interaction is given by 


Q. x = aQxH[-{l + AYy,{V,-B)Y^. (18) 


From (10) and (18), we obtain that the perturbed inter¬ 
action can be written as 


SQx = 6F[Qx{a) + HQ'^{a)5a , (19) 

while for the momentum transfer potential we hnd that 
Fx = Qx{Vc-V) = -Fc provided Q^\\u^, avoiding in 
this way the violation of weak equivalence principle [28j . 
Notice that SQx contains two different kinds of terms, 
namely SQx receives contributions from the spatial vari¬ 
ation of H and the spatial variation of the scale fac¬ 
tor. In order to carry on and explore the impact of this 
new interaction on the parameter space, it is essential 
to modify or include several routines in the numerical 
code called CAMB. Given that this code is written in 
the synchronous gauge we need to handle all these ex¬ 
pressions into this gauge characterized hy A = B — 0. 
Replacing both conditions into gauge transformation for 
metric variables indicates that scalar perturbations can 
be treated in terms of two functions: rjT = — Hl 

and h]^ = [1H]- From the later fact, with the help of 

the metric ([^ , we calculate SF[ and Sa in the synchronous 
gauge as 


[5a]syn = -a{T]T - T-La) , 

[57P]syn = -^[?)T + a(H2-H)], (20) 


where cr = {hi, -I- 377t)/ 2A:^ and the overdot refers to con¬ 
formal time derivative. Combining Eqs. (20) and Eq. 
(19) in the synchronous gauge, we have 


SQx = ^[PT + P'HvT + cT{n^ - -n)]. ( 21 ) 


To understand the cosmological constraints performed 
with redshift space distortion measurements we must give 
further insights about the dynamical evolution of dark 
matter perturbations. We need to examine if the ex¬ 
change of energy-momentum in the dark sector may al¬ 
ter the standard Euler and continuity equations for dark 
matter, leading to a detectable signature about its dy¬ 
namic in relation with the behavior of galaxies. Our 
starting point is to calculate the evolution equations for 
the perturbation of cold dark matter in the synchronous 
gauge |2H]: 


L + + ^ = -[SQc-ScQcy (22) 

2 Pc 

Oc + Ocn = -[Qc{o-ec)-k^Fcy ( 23 ) 

Pc 
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where we used that the adiabatic sound speed defined 
as = vdPc = Wc + WcI{PcIPc) vanishes, Wc = 0, 
and the physical sound speed in the rest frame, namely 
cL = {Spc/Spc)rest frame = 0 for cold dark matter, and 
the density contrast is defined as Sc = Pc^Pc/Pc- For the 
definition of 9c and 9 see our previous article [28] . 

After decoupling, the baryons do not interact with the 
dark sector so its Euler and continuity equations within 
the synchronous gauge are given by 


4 


9b + — — 0 , 

9b + 9bn = 0. 


(24) 

(25) 


In our case, we find [Qc(9 — 9c) — = 0 provided 

the momentum transfer potential is —dFc = Qc{9c — 9). 
This implies that the Euler equation for dark matter is 
not altered in relation with the non-interacting case. As a 
result , we ob tain that the difference between Euler equa¬ 
tions ( 23p5 ) leads to {9c — d) = —T~L{9c — 9b), and there¬ 
fore we can assure that there is no velocity bias within 
this model [41]. This fact pinpoints that dark matter ve¬ 
locity is not directly affected by the interaction, which 
means that dark matter follows geodesic and there is not 
violation of the wea k eq uivalence principle |42j . 

Combining Eq. (22) for d along with the Einstein 
equations in the synchronous gauge |40|, we arrived at 
the second-order evolution equation for dark matter per¬ 
turbation: 


Sc + US, 




Pc'H/ 


47rGa^ 


SbPh T ScPc ^ 


1 - 


{Pc ^Qc)- 

AirGapc 




The growth of dark matter perturbations are affected in 
a number of ways. Firstly, the modification of the dark 
matter contrast density enters through the expansion his¬ 
tory encoded in the Hubble expansion and cold dark mat¬ 
ter density, which will lead to a faster or slower cluster¬ 
ing of matter particles according to the values taken by 
d and /3. Secondly, the friction term of Eq. (26) is al¬ 


tered by the extra contribution, aQdPcT~L, implying that 
the rate growth will depend both on Qc and the back¬ 
ground evolution described by %. Thirdly, the cold dark 
matter particle feel and additional force that is mediated 
by dark energy particle described by an effective gravi¬ 
tational constant. As cold dark matter particles have a 
mass which varies with time, the uncoupled baryons will 
feel this effect through a changing gravitational poten¬ 
tial. It is clear that the extra-term in Eq. (26) is related 
with the fact that SQc cannot be recast as SpcQc pro¬ 
vided our interaction is no proportional to Pc- At this 
point, we must highlight the importance of dark mat¬ 
ter perturbations and baryon perturbations, both con¬ 
tributions are essentials for understanding the formation 
of the large-scale structure of the Universe. More pre¬ 
cisely, the growth data involve the linear perturbation 


z 


survey 

0.067 0.42 ± 0.06 

6dFGRS (2012) 

0.17 

0.51 ±0.06 

2dFGRS (2004) 

0.22 

0.42 ±0.07 

WiggleZ (2011) 

0.25 

0.39 ±0.05 

SDSS LRG (2011) 

0.37 

0.43 ±0.04 

SDSS LRG (2011) 

0.41 

0.45 ±0.04 

WiggleZ (2011) 

0.57 

0.43 ±0.03 

BOSS CMASS (2012) 

0.60 

0.43 ±0.04 

WiggleZ (2011) 

0.78 

0.38 ±0.04 

WiggleZ (2011) 

0.80 

0.47 ±0.08 

VIPERS (2013) 


TABLE I: Compilation of fas{z) data points obtained from 
several galaxy surveys using RSD method. 


of the overall growth factor dm in terms of the function 
fm = din dni/dln a along with the r.m.s density contrast 
within an sphere of radius i?8 = 8Mpc h“^, which is re¬ 
lated with the matter power spectrum. As explained be¬ 
fore, the measurements of redshift space distortions lead 
to a bias-independent quantity called f{z)(7s{z) [H] and 
is used it as a powerful tool for improving the cosmolog¬ 
ical constraints beyond the background probes based on 
supernovae data. It is possible to use these data in our 
cosmic constraint but we must first modify the CAMB 
code for extracting the theoretical value of f{z)as{z), and 
such procedure involves to re-write different routines. On 
one side, (Ts{z) is given by 


crsiz) 


1 

2tt'^ 



dkdwi{k)Pik, 



(27) 


where Ws{k) is the Eourier transform of the top-hat win¬ 
dow function, with an interval defined by Rs, and P{k, z) 
refers to the matter power spectrum. It can be seen that 
in our model the overall growth matter rate is 


^ PcS'c + PbS'b dPcpa^ 

R-iSc Pb) Pc Pb \ / 


(28) 


where the interaction induces an extra term which plays 
the role of a relative weight between cold dark matter 
contrast density and total matter contrast density. More¬ 
over, ([28]) shows that the interaction reduces the overall 


growth matter rate for positive coupling constant, d > 0- 
This difference therefore must leave to a detectable signa¬ 
ture, which allows differentiation between previous works 


III. CONSTRAINT AND RESULTS 

We perform a statistical estimation of the cosmic 
parameters by using the Markov Chain Monte-Carlo 
method with help of public code CosmoMC |36| and 
CAMB code |37|. In order to do so, we add a new module 
to CosmoMC package for taking into account the value 
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FIG. 1: 68 % and 95 % constraints on w^-, and /3 from the combined analysis made with the Planck 2013+WMAP9+JLA+BAO+RSD+HST 
data. The ID marginalized posterior distribution of w^, and /3 are also shown. 
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FIG. 2: Typical effects of the coupling constant on the theoretical CMB temperature power spectrum. We exhibit versus 
multipole for different values of the coupling constant. 
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FIG. 3: Typical effects of the j3 exponent on the theoretical CMB temperature power spectrum. We show versus multipole 
for different power-law interactions. 
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FIG. 4: The total matter power spectrum at 2 = 0 as the coupling constant varies from 0 to 0.9. Baryon acoustic oscillations generate wiggles 
the matter power spectrum. 
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FIG. 5: The total matter power spectrum at 2 = 0 as the 0 exponent varies from 0 to -1. Baryon acoustic oscillations generate wiggles in the 
matter power spectrum. 
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Parameters 

Mean with errors 

Best fit 


0.02215594+”;“fg 

0.02222252 


0.118649101E]:“« 

0.11845120 

1006MC 

1 fi/l 1/I 

< -0.00057134 

1.04152800 

T 

0.08910617l°°l"3l®^^ 

0.08779789 

Wx 

1 1 nQ'7/i on Q”l"0.08205012 
— i.iZO/06244387 

-1.14908200 


U.044U04Z ( -0.18907353 

0.32413630 

P 

U.OUOOOZOZ_Q 40^23857 

-1.06008500 

ln(10i°A,) 

3.08464588l°;g2425M5 

3.08247800 

Us 

0.96269367l°“°®24026 

0.96447760 

Hq 

63.68287359l^j««4i6 

65.75439000 


0.64884118l[;;°g^2i75 

0.67314830 


n 1 c:qo9+0-02299549 

U.0G11000Z_q Q4222199 

0.32685170 


0.14145019l[]°[]i«I?^t 

0.14131890 

ers 

0.769852121°;“?“®“ 

0.78026930 

Age/Gyr 

1 77oqi::p:i::9+0.03828760 

LO.i ' Q3gQ4Qg0 

13.78876000 


1090.0590572112;?“®““ 

1089.95200000 

r* 

144.94953514l°;??“26« 

144.94740000 

0. 

1.04165139l2;22256?87 

1.04174300 


TABLE II: Statistical results from the 

global fitting performed with the Planck 
2013+WMAP9+JLA+BAO+RSD+HST data. 


of fm obtained with the CAMB code, in this way, we cal¬ 
culate the theoretical value of f{z)as{z) at different red- 
shift, compare with their observational value, and build 
the likelihood for performing a statistical estimation 
on the parameter space [43]. Two comments are in order 
concerning the PPF “parameters”. When cp is taken as 
free model parameter one finds that its posterior prob¬ 
ability distribution has a flat shape, showing that any 
value between zero and the unity are physically admissi¬ 
ble |2H|, in fact, a similar analysis can be obtained for the 
function [3T|. The list of the cosmological data used 
for our parameter estimations (cf. |5S] for further details 
about the methodology) are: 

• SNela, JLA data: 740 supernovae samples from low 
redshift z = 0.02 to large one z ~ 1, obtained from 
the joint analysis of SDSS II and SNLS [3|. 

• CMB, WMAP9 -I- Planck: multipole measurements 
obtained by WMAP9 team |2] and Planck satellite 
SI, [5], 0. WMAP9 project involves the measure¬ 
ments of Atacama Cosmology Telescope (ACT) at 
high multipoles, £ S [500,10000], along with the 
South Pole Telescope (SPT) observations which re¬ 
ported data over the range i G [600,3000]. Planck 
survey performed measurements over a complemen¬ 
tary zone, i G [2,2500]. 

• BAG, DR9 - DR7 - 6dFGS: the 6dFGS mis¬ 
sion reported (0.106) [3, SDSS-DR7 measured 
dz(0.35) |8], SDSS-DR9 exploration led to dz(0.57) 
[Sj, and diverse measurements from the Wiggle Z 
dark energy survey reported ^^(O.dd), (iz(0.60), and 
4(0.73). 

• HST, Hubble: Gaussian prior for Hq. 


• RSD, Growth data: measurements of the quantity 
f{z)as{z) at different redshifts unify the cosmic 
growth rate / and the matter power spectrum erg 
normalized with the co-moving scale 8h“^ Mpe, in 
a single quantity which includes the latest results of 
galaxy surveys such as 6dFGS, BOSS, LRG, Wiggle 
Z, and VIPERS [see Table ([!])]. 

The parameter space is given by [55] 

V = Vtch?, 1006»mc, ns, ln(10^°As), r, /?), 

(29) 

and each data base is taken as independent one, so the 
likelihood distribution is 


Xtotal — XSNe + XbAO + XCMB + XhST + XrSD- (30) 


Let US analyze the outcome of our sta¬ 
tistical analysis obtained by combining all 
Planck2013 -b WMAP9 -b JLA -b BAG -b HST -b RSD 
data [see Table 0 ]- The equation of state of dark energy 
is given by Wx = — 1.12874223lQ;Qg2443g7 at Icr level, 
such values confirm that the observational data prefer 
a phantom-like equation of state over a non-phantom 
one [see Fig. 0]. Regarding the amount of dark matter 
at present time we find that the best fit to all data 
combined: £lrn = 0.35115882to;o4222i99 level, 

while, the Planck -b WP data leads to flm = 0.315to;oi8: 
then the disagreement is not bigger than 0.11%, 
moreover, comparing with the Planck 2015 data we 
obtain a relative difference near 0.22% for the Planck 
EE-blowP [JS]. The best fit for dark energy amount 
is fix = 0.64884118^Q;g2299489 level, showing a 

deviation of 0.05% in relation with the value reported by 
Planck mission [^. This result is utterly consistent with 
the interacting cosmology explored here because the 
transfer of energy goes from dark energy to dark matter, 
then we can expect that such mechanism increases the 
current value of dark matter. Further, the dynamical 
probe included in the statistical analysis clearly shows 
that the overall growth rate of matter (281 is quite 
sensitive to the interaction proportional to the current 
dark matter amount. In this regard, the joint analysis 
of Planck20I3 -b WMAP9 -b JLA -b BAG -b HST -b RSD 


data [see Table 0 and Fig. 0 ] shows that a large cou¬ 
pling constant cannot be ruled out in our model provided 
= 0-34403427lo;^g9O7353 at I ct level. We stress that 
this result depends on the addition of growth rate data 
and the interaction itself. More precisely, if one considers 
only geometric probes for performing the statistical test 
with CMB -b BAG -b SNe data, the coupling constant 
remains lower, near 0(10“^) m- Despite the inclusion 
of the growth rate data gives a non zero deviation 
from the previous result, the strength of the interaction 
still remains small, that is, = O-OOWO^aoooso 
level [28]. Notice that the inclusion of a perturbed 
expansion Hubble contribution in the standard pertur¬ 
bation theory does not prefer large interaction even 
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if growth rate data plus geometric probes are taken 
into account [44) . Therefore, the key point for under¬ 
standing the large value of coupling constant reported 
here is linked with some cooperative effects coming 
from the inclusion of perturbed expansion Hubble rate 
along with spatial variation of scale factor within the 
interacting PPF method. Besides, the combined data 
of Planck2013 -b WMAP9 -f JLA -b BAO -b HST -b RSD 
prefer a power law form with negative exponent, namely 
P = —0.50863232+Q 4 Qg 23 g 5 Y at 1 cr level, indicating 
that the interaction is stronger in the early universe 
and weakens at present time. However, notice that 
our finding does not mean that the interaction has a 
vanishing contribution in the latter epoch. 

The impact of the interaction can be traced by inspect¬ 
ing the changes introduced in the cosmic background ra¬ 
diation and clustered matter power spectrum, thus we ex¬ 
plore how these quantities are affected by the two param¬ 
eters which describe the interaction itself: and /3. Fig. 
([^ shows the CMB temperature anisotropies in terms of 
the multipoles for several values of the coupling constant 
with G [0,0.9]. Firstly, the position of the acoustic 
peaks shifts to right by increasing the coupling constant. 
The peaks are located at where the acous¬ 

tic scale is defined as 6a = s{zdec)l‘rzj,^^ with r known as 
the comoving angular diameter distance, Zdec — 1090 in¬ 
dicates the epoch when baryons decouple from photons, 
and the sound horizon is given by 


s{z) 




Cs(a) 


da 

a^H{a)' 


(31) 


where th e sound s peed for the photon-baryon fluid is 
Cs = cj -\/3(l -b i?) and R = Since the cou¬ 

pling constant is positive, the exchange of energy goes 
from dark energy to dark matter, increases linearly the 
amount of dark matter at early time, affecting the sound 
horizon at the end, which modifies the value of the first 
peak located at £i ~ 7r/6*A as it can be seen from Fig. 
([^ by comparing the red-line (^c = 0.34) with the black¬ 
line (^c = 0). A secondary effect refers to the amplitude 
of CMB power spectrum at lower multipoles, it turns 
out that increasing the coupling constant shows a de¬ 
viation from the vanilla model, so we believe that such 
phenomenon can be properly explored in the future by in¬ 
cluding the galaxy-ISW cross-correlated power spectrum 
m- Besides, Fig. ([^ shows a shift to right in the acous¬ 
tic position peaks as one moves from a vanishing expo¬ 
nent to a negative exponent equal to —1 while the other 
parameters are fixed to our best cosmology. This effect 
can be understood easily by taking into account that the 
extra-term in the cold dark matter density increases from 
a constant value to a~^, then it amplifies the amount of 
dark matter at early times. It is clear that changes in 
the background expansion history of the Universe and 
modifications to the growth rate of matter perturbations 
will affect the matter power spectrum. In Fig. we 
plot the matter power spectrum at z = 0 generated us¬ 
ing CAMB for different values of the coupling constant. 


while. Fig. ([^ depicts the behavior of P{k) when the 
exponent P changes from zero to —1. Fixing all param¬ 
eters to our best fit cosmology except or the absolute 
value of p will have almost the same kind of impact in 
the matter power spectrum. The growth of density per¬ 
turbations are significantly affected by the variation of 
P or This is because the amount of dark matter is 
bigger than the fraction of dark energy at early times 
when > 0 and /3 < 0. As a result the epoch of matter- 
radiation equivalence happens earlier than the standard 
case so only the very small scale perturbations on scale 
k > keq have time to enter the horizon and grow during 
the radiation domination. Hence, the most obvious effect 
is that the turnover in the power spectrum (which is as¬ 
sociated with the scale that entered the horizon when the 
Universe became matter dominated) is placed on smaller 
scales. In addition, the amplitude of P{k) is considerably 
affected as can be noticed from Figs. We must 

stress that baryons also leave their impact on the matter 
power spectrum due to their coupling with the photons 
before recombination. Indeed, they produce a series of 
wiggles in the power spectrum [see Figs. (§] , which 
are related with the existence of a well defined peak in 
the correlation function of galaxies, ^(r), placed at 150 
Mpc [1]. 


IV. SUMMARY 


To avoid large-scale instability at early times associ¬ 
ated with the growth of dark energy perturbations, we 
have employed the interacting PPF method for exam¬ 
ining a new kind of interaction which includes spatial 
variation of Hubble function along with spatial variation 
of scale factor. Interestingly enough, we found that the 
transfer of energy-momentum can flow from dark energy 
to dark matter only, allowing to preserve dark energy and 
dark matter densities as positive quantities for all times. 
Concerning matter perturbations, we have shown that 
the overall growth rate of matter is strongly affected in 
three different ways. Firstly, the dark matter density and 
the Hubble rate at background level both deviate from 
the concordance model. Secondly, an extra-term related 
with friction in (26) reveals that the local Hubble func¬ 
tion is modified, and thirdly the gravitational constant 
is amended also. The aforesaid effects are essentials for 
understanding the cosmological constraint based on the 
overall growth rate function fm = din dm/dln a provided 
is connected with the measurements of redshift space 
distortions through a bias-independent quantity called 

f{z)Mz) [E]- 

Regarding the observational side, we per¬ 
formed a global fitting based on geomet¬ 
ric and dynamical probes by combining 
Planck2013 -b WMAP9 -b JLA -b BAO -b RSD -b HST 


data [see Table (|TT])]. The joint analysis showed that 
a phantom equation of state is preferred at 1 u level 
and the amount of dark matter differs from the Planck 
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2013 estimation in 0.11% while the Planck 2015 doubles 
this difference after the Planck EE+lowP data is used 
into the analysis [IS]. Such disagreement can be easily 
understood by taking into account that the interacting 
mechanism described here only allows the exchange 
of energy from dark energy to dark matter. Another 
important ingredient concerns to the particular kind of 
power law form selected by the data and the strength 
of coupling constant. More precisely, the joint analysis 
of Planck2013 + WMAP9 + JLA + BAO + RSD + HST 
data [see Table 0 and Fig. 0 ] showed that a large 
coupling constant cannot be ruled out in our model 
provided, = 0-34403427lo;i[g9g7353 at 1 cr level. At 
this point, we must note that previous estimations 
give a coupling constant near 0(10“^) [37], [35], even 
though the perturbed expansion Hubble rate term is 
included into the analysis [44]. All the combined data 
(Planck2013 + WMAP9 + JLA + BAO + RSD + HST) 
singles out a power law form with negative exponent, 
P = —0.50863232to;4og23857 at 1 O' level, pinpointing the 
crucial effect of the interaction at early times. 

We have examined how the interaction affects the tem¬ 
perature CMB spectrum profile, focusing in the position 
of acoustic peaks and their amplitudes. A similar analy¬ 
sis was carried out for the matter power spectrum where 
the target was to explore the position of turnover scale. 
Varying only the conpling constant, we found that the 
position of the first peak in the TT power spectrnm is 


shifted because the fraction of dark matter increases lin¬ 
early with [ see Fig ([^]. As the /3-exponent varies only 
a similar effect can be pinpointed [see Fig ([^]. Besides, 
the amplitnde of the TT CMB spectrum at low multiples 
is sensitive to the value taken by the coupling constant, 
indicating that the integrated Sachs-Wolfe must be in¬ 
cluded in a near future analysis m- We also found an 
interesting effect in the matter power spectrum by vary¬ 
ing or p. It turned out that the scale of matter and 
radiation equality is shifted further from the present era 
when the amount of dark matter is increased by varying 
or P, and therefore a shift of the turnover point toward 
to smaller scale is detected [see Figs. Q- ([^]. 
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